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A ROTATING MOLECULAR DISK TOWARD IRAS 18162-2048, 
THE EXCITING SOURCE OF HH 80-81 

M. Fernandez-Lopez^'', J.M. Girart^, S. Curiel\ Y. Gomez'', P.T.P. Ho^"', and N. Patel"' 



O 



< 

Oh- 
I 

o 

(N 
> 

m 

o 



% 



ABSTRACT 

We present several molecular line emission arcsec and 
subarcsec observations obtained with the Submilhmeter 
Array (SMA) in the direction ol the massive protostar 
IRAS 18162-2048, the exciting source of HH 80-81. 

The data clearly indicates the presence of a compact 
(radiusw 425-850 AU) SO2 structure, enveloping the 
more compact (radius<150 AU) 1.4 millimeter dust emis- 
sion (reported in a previous paper). The emission spa- 
tially coincides with the position of the prominent ther- 
mal radio jet which terminates at the HH 80-81 and 
HH 80N Herbig-Haro objects. Furthermore, the molecu- 
lar emission is elongated in the direction perpendicular to 
the axis of the thermal radio jet, suggesting a disk-like 
structure. We derive a total dynamic mass (disk-like 
structure and protostar) of 11-15 Mq. The SO2 spec- 
tral line data also allow us to constrain the structure 
temperature between 120-160 K and the volume den- 
sity > 2 X 10^ cm~'^. We also find that such a rotating 
flattened system could be unstable due to gravitational 
disturbances. 

The data from C^^O line emission show a dense core 
within this star-forming region. Additionally, the H2CO 
and the SO emissions appear clumpy and trace the disk- 
like structure, a possible interaction between a molecular 
core and the outflows, and in part, the cavity walls ex- 
cavated by the thermal radio jet. 

1. INTRODUCTION 

The powerful radiation from high-mass protostars 
was the main obstacle for understanding the forma- 
tion o f mass i ve stars (see e.g. iLarsoii fc Starrfieldlll971l: 
Kall^ [T97I lYorke fc Kruegell Il977t IBeech fc Mitalas 
19941 : see also the recent review of iZinnecker fc Yorke 
2007| ). It was thought that the strong radiation 
pressure from the protostar would stop the accre- 
tion, hence preventing the formation of stars with 
masses over '^10 Mq. Alth ough this theoretical prob- 
lem can be overcome (e.g., iWolfire fc Cassinellil 119871 : 
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Fig. 1. — Contour image of the zero order moment superposed 
on the first order moment colour image of the C^^O 2—1 line emis- 
sion (combined compact and extended configuration data with a 
Gaussian uv taper) toward IRAS 18162-2048. Contour levels are 
-9, -4, 4, 9, 15, 30, 45, 60, 75 and 90 times 0.10 Jy beam"! km s-\ 
the rms of the image. The units of the scale bar are in km s~^ . 
The molecular emission is integrated over velocities between 9 and 
18 km s~^. The triangles mark the position of the millimeter 
sources MMl and MM2, and the position of the molecular core, 
MC. The dashed arrows show the directions of the NW and SE 
outflows which seem to originate at the MM2 position and also the 
NE-SW thermal radio jet, launched from MMl. The synthesized 
beamsize is shown at the bottom right corner and the primary 
beam is indicated by the dashed circle. 
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Krumholz et al.ll2005[ ). it remains yet unclear how these 
stars are actually formed. There are several hypothe- 
ses for the formation process. Perhaps the most im- 
portant ones are: (1) large accretion rates (3 or 4 or- 
ders of magnitude greater than those observed in low- 
mass protostars) through circumstellar disks (|W almslc3 
19951 : 1 Jiiina fc Adams! [1995 IMcKee fc Tanl l2003i: .Zhai^ 
20051 iBaneriee fc Pudritzl l2007t) . (2) competitive ac- 



cretion between small protostars of the same cluster 
that could result in eventual mergers U 3onncll ct al 
1991 iClarke fc Bonncll 200l iMoeckel fc Clarkd l201ll 
Baumgardt fc Klesscn .20111) . and (3) accretion of ion- 



ized material, even after the formation o f a compact 
HH region generat ed by the protostar (jKetol 120021 : 
IKeto fc Woo"dl[200l . 

In recent years, several works have provided evi- 
dence for the pr esence of coUimated jets (HH 80-81, 
Marti et al.l[T9M Cep A. ICuriel et aU&m. IRAS 23139 
Trinidad et all [2001 IRAS 16547. [Rodriguez et al.l 



20081). molecular outflows (see a summarv in lZhanel[2005[ : 
Arce et al1[2007D . and flattened structures of dust and 



gas (some with rotation and some with infall signposts; 
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Spectral observations 



-20 47'40" 



Transition 


Inline 


Eu 


Data 


Synthesized 


Beam 
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SO2 163,13-162,14 


214.689380 


147.9 


LR 


8'.' 3 X 


3'.'1 


34.2 


0.065 










ROB 


1'.'3 X 


o'.'g 


59 


0.085 










HR 


0'.'7 X 


0'.'4 


13.2 


0.085 




SO2 176,12-185,13 


214.728330 


229.1 


LR 


8'.'3 X 


3'.'1 


34.2 


0.065 






SO2 64,2-73,5 


223.883569 


58.6 


LR 


7'.' 9 X 


3'.'0 


34.2 


0.065 






SO2 202,18-193,17 


224.264811 


207.9 


LR 


7'.' 9 X 


3'.'0 


34.2 


0.065 






SO2 132,12-13l,l3 


225.153702 


93.1 


LR 


7'.' 9 X 


2'.' 9 


34.2 


0.065 
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1'.'2 X 


O'.'O 


58 
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O'.'T X 


0'.'4 


12.9 


0.080 




SO 78-77 


214.357004 


81.3 
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8'.' 3 X 


3'.'1 


34.2 


0.065 




SO 55-44 


215.220653 


44.1 
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3'.'1 


34.2 
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Ci^O 2-1 


224.714385 


16.2 


LR 


7'.' 9 X 


2'.' 9 


34.2 


0.075 


(a) 








TAP 


4'.'0 X 


2'.'1 


37.8 


0.080 


ID 


H2CO 3i,2-2i,i 


225.697775 


33.5 


LR 


7'.'9 X 


2'.'9 


34.2 


0.090 


101 








ROB 


1'.'2 X 


O'.'g 


58 


0.100 


El 








TAP 


4'.'0 X 


2'.'1 


37.2 


0.095 |fi|5|4| 



Note. — Data extracted from JPL catalogue. The fourth column Data is a code giving the 
resolution of each map. LR means low angular resolution data, HR, high angular resolution 
data, ROB, combined images with the robust set to 0.3 and TAP, combined images with a 
taper applied and restricted in the uvplane (see text, section §2). 

(a) The compact configuration data of the C^^O 2-1 is not shown in any figure. However, it 
is used to obtain the velocity, density and mass of the dense core of IRAS 18162-2048 (see 
section §3.1). 
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Fig. 2. — Contour image of the zero order moment of the SO 55- 
44 line emission (combined compact and extended configuration 
data with a Gaussian uv taper) superimposed on a K-band near 
infrared image of the refiection nebula in the IRAS 18162-2048 re- 
gion (courtesy of Thomas Geballe; seeJAspin & Geballc 1992 ). The 
infrared polarimctric study carried out bv iYarnashita et anill987l ) 
demonstrated that the scattered light from the infrared refiection 
nebula occurs on grains at the walls of a parabolic cavity rather 
than on grains inside the whole outfiow lobe. Contour levels are -9, 
-4, 4, 9, 15, 30, 45, 60, 75 and 90 times 0.17 Jy beam-^ km s-\ 
the rms of the image. Symbols are as in Fig. [l] 



e.g., Cep A, IPatel et all [200l IRAS 201 26. ICesaronil 
20Cia AFGL 2591 Ivan der Tak et al.ll2006t W51 North 
Zapata et al.ll2008l: IR AS 165 47,'Franco-Hcrnandc z et alj 
20091; W33A-MM1, iGalvan-Madrid et al. 2010), sur- 
rounding high-mass protostars. The evidence leads to 
the interpretation that massive star formation is analo- 
gous to low-mass star formation, that is, via accretion 
from a flat rotating disk, with a jet of ionized material, 
and with an associated molecular outflow. However, the 



physical characteristics of the possible accretion disks 
and their associated outflows have not been well char- 
acterized yet. Hence, the formation process of massive 
stars remains an open issue. There are several prob- 
lems that hamper the analysis of massive star formation 
regions (MSFRs). First, massive protostars form in clus- 
ters, which together with projection effects, make these 
regions difficult to interpret. Moreover, the MSFRs are 
found at larger distances than low-mass star formation 
regions, with typical distances of 2-5 kpc. Thus, to get 
an insight of the closest regions to the protostars, where 
accretion disks are expected, the highest angular resolu- 
tion of the current telescopes is needed. Even so, most of 
the dust and gas structures detected around high-mass 
protostars are ~ 1000 AU. This is the smallest diameter 
that an interferometer can resolve with a 1" angular reso- 
lution at those distances. At these scales, such structures 
could harbour systems of several protostars, that could 
be forming high mass stars by other mechanisms such as 
mergers. In fact, the nearest massive protostars (Orion 
BN/KL and Cepheus A HW2), show complex scenar- 
ios, with a possible merger of protostars (Orion BN/KL, 
iZapata et al.l l20()9'). or close passages between protostars 
(Cep A HW2]Cunningham et al. 200i). Therefore, with 
the current interferometers, only a few MSFRs may be 
studied with adequately high angular resolution. 

To demonstrate the existence of a disk around a 
massive protostar, it is not enough to detect its 
(sub)millimetcr dust continuum emission; its kinemat- 
ics must also be measured via the molecular line emis- 
sion associated with the gas-phase chemistry developed 
after the evaporation of the molecules from the ice man- 
tles (JHatcheiriFiD [Till [CeiaromiriD [200^ How- 
ever, searching for molecular tracers of disks in MS- 
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Fig. 3. — Channel velocity image of the SO 65-44 extended emission (combined compact and extended configuration data with a Gaussian 
uv taper) from the central region of IRAS 18162-2048. The emission is averaged in velocity bins of 1.14 km s~^. The velocity of each 
channel is indicated in the top left corner of each panel in km s~^ units. The cloud velocity is at about 11.8 km s~^. The contour levels 
are -6, -4, 4, 6, 9, 13, 18, 24, 31, 39, 48 and 58 times 0.08 mjy beam" , the rms of the image. Symbols are as in Fig. [T] The synthesized 
beam size is shown in the bottom left-hand corner of the most redshifted channel. The blueshiftod emission appears to trace mainly the 
walls of the cavity excavated by the radio jet from MMl. At the position of MMl the image shows emission probably coming from the 
rotating disk/ring, absorbed in the bluest channels by the extended emission. The redshifted channels show also compact emission toward 
the position of MC. 
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Fig. 4.— SO 55-44 (top panel) and H2CO 3i,2-2i,i (bottom 
panel) spectra toward the center position of MC (black line). Spec- 
tra toward the clump located to the north-west of MC (green line) 
are included to compare with that of MC. The dotted line marks 
the systemic velocity (11.8 km s~^). The intensity is in Kelvin. 
The spectra are extracted from the combined compact and ex- 
tended configuration data with a Gaussian uv taper. 



FRs is a complex task. There are molecules that 
show emission from the envelope and the disk simul- 
taneously. Other molecules show optically thick emis- 
sion, thus complicating the kinematic study of the 
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Fig. 5. — Contour image of the zero order moment of the 
H2CO 3i_2-2i,i line emission (combined compact and extended 
configuration data with a Gaussian uv taper) overlaid with a K- 
band near infrared imago of the IRAS 18162—2048 region (as shown 
in Fig. [2J. Contour levels are -9, -4, 4, 9, 15, 30, 45, 60, 75 and 90 
times 0.33 Jy beam^^ km s^^, the rms of the image. Symbols are 
as in Fig. [Tl 

disk. On the other hand, S-bearing species (such 
as H2S, SO, SO2, CS, OCS...) could be intimately 
linked with the evaporation process of the disk surface 
(lCharnlevlll997t iHatchell et all 1199 81: Ivan der Tak et al.l 
I2003t iMartin-Pintado et all |2005[) . becoming good trac- 
ers of the dynamics of the innermost parts of the high 
mass protostars. Recently, several papers have been pub- 
lished on the detection of S-bcaring s pecies in disks and 
other warm gas-structures of MSFRs (Ivan der Tak et alj 

120091 : 



20061 iJimenez-Serra et alJ l2007t iKlaassen et al 



Franco-Hernandez et al.ri2009t iZapata et al.l 120091 ) . In 
particular, SO2 transitions, ubiquitous within the 
(sub)millimetcr range, show a very compact nature, sug- 
gesting their close association with circumstcUar struc- 
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Fig. 6. — Channel velocity image of the H2CO 3i,2-2i,i extended emission (combined compact and extended configuration data with 
a Gaussian uv taper). The emission is averaged in velocity bins of 1.14 km s^^. The velocity of each channel is indicated at the top left 
corner of each panel in km s~^ units. The cloud velocity is at about 11.8 km s~^. The contour levels are -6, -4, 4, 6, 9, 13, 18, 24, 31, 39, 
48 and 58 times 0.095 mjy beam"^, the rms of the image. Symbols are as in Fig. [T] 



tures. 

The GGD27 complex is an active star forming re- 
gion located at a distance of 1.7 kpc. It shows a 
spectacular (5 .3 pc long) and highly coUimated ther- 
mal radio jet (jMarti et al.lll993l 119951 I1998D. embedded 
within a powerful CO bipolar outflow ( Yamashita et al] 
Il989t iBenedettini et all I2004D The jet has a position 
angle of ^ 21° (|Marti et al.l [19931 iMarti et al.l [TMl) 
and ends at two southern and very bright Herbig- 
Haro objects (HH 80 and 81, originally discovered by 
iReipurth fc Grahainlll988D and at a rad io source to the 
north (HH 80 North. iGirart et al.l [1991) . Linear polar- 
ization has been detected in the radio emission from this 
jet, indicating the presence of a magnetic fie ld coming 
from the disk (jCarrasco-Gonzalez et al.ll2010f ).The cen- 
tral part of the radio jet has a bright far-infrared coun- 
terpart (IRAS 18162-2048), that implies the presence of 
a luminous young star (~ 2 x 10^ \jc.-^ or a clu ster of stars 
(jAspin & Geballel [T992: Stc cklum etall [19971 ) . 

Submillimeter and millimeter wavelength observations 
of the central part of the thermal radio jet have re- 
vealed two du sty sources, MMl and MM2, sep arated 
by about 7" (|G6mez et al.l [200l iQiu fc Zhanil[2009[ ). 
which arc apparently in very different evolutionary stages 
(|Fernandcz-L6pez et al. 20111 hereafter Paper I). MMl, 
the south-western source, coincides with the origin of the 
thermal radio jet and has a dust temperature of 109 K. 
It has a weak extended envelope, possibly surrounding a 
very compact (R<150 AU) disk-protostar system. The 
mass of the disk derived from the continuum emission is 
~ 4.1 Mq. Such a massive disk could have an extremely 
high accretion rate (10""^ — 10~^ Mq yr~^). However, 
the case for an accretion disk in MMl has not been un- 
ambiguously confirmed yet. The bolometric luminosity 
obtained from fitting the SED of high angular resolu- 
tion data is > 3300 Lq, similar to that of a Bl Zero 
Age Main Sequence (ZAMS) star. However, the mas- 
sive disk of MMl (compared to those of low-mass pro- 



tostars), the presence of the powerful outflow and the 
youth of the protostar (it has not developed a com- 
pact H II region yet), indicate that MMl could become 
a BO or a more massive star. MMl is also associated 
with compac t emission from several high d ensity tracers, 
such as CS ([Yamashita et al.l Il99i[ ). SO (jGomez et al.l 
[200l. CH. OH. CH3CN, Hi^CO. OCS, HNCO and SO2 
(|Oiu fc Zhang 2009) . 

The physical characteristics of M M2 are typical of th e 
Class low-mass protostars (see e.g.. lAndre et al.iri993() . 
except for its total mass of at least 11 Mq (4-5 times 
mor e mass than a typical low-mass Class source, see 
e.g. IGirart et all [2006: Rao ct al. 200§). MM2 appears 
as an extended source when observed with low angular 
resolution at 1.4 mm, but it is separated into two com- 
pact sources when observed with high angular resolution 
(Paper I). The stronger component in MM2 spatially co- 
incides with very weak free- free emission, suggesting that 
it could be a high/intermediate mass protostar. The 
weaker component in MM2 could be a pre-protostellar 
core. Previous CO observations (|Qiu fc Zhang] |2009() 
have associated the position of MM2 with the origin of a 
young outflow running to the east. In addition, another 
possible CO outflow, running north-west of MM2 could 
also have originated close to this position. 

There is also a po ssible molecular core (hereafter MC, 
lOiu fc Zhang|[2009f ) located about 4" to the north-west 
of MM2. Although MC is traced by several molecules 
(CH3CN, H13C0, OCS y HNCO. lOiu fc Zhand [20091 
and a Class I CH3OH maser. [Kurtz et al.ll2004f )rit is not 
associated with a bright and compact millimeter source 
(Paper I). Th e diverse line emissio n from MC has been 
interpreted by lOiu fc Zhang] (|2009() as the emission from 
a hot core warmed by the radiation of a young and mas- 
sive protostar. 

Here we present an analysis of the molecular emission 
of several lines (mostly S-bearing species) detected in 
the central region of IRAS 18162-2048 using the Submil- 



limcter ArrajQ (SMA: lHolFaL|[200l . The SMA allowed 
us to detect several transitions simultaneously, some of 
which seem to trace a circumstellar disk/ring around 
MMl, with subarcsecond resolution (~ 0'.'5), permitting 
a beam size of ^ 850 AU at the source distance. This 
resolution is necessary in order to perform a reliable kine- 
matical analysis of the motions of the circumstellar disks 
at distances greater than 1-2 kpc. 

The continuum emission of the observations was pre- 
sented in a recent publication (Paper I). In section 2 we 
briefly describe the observations, while in section 3 the 
results of the observations are given. Section 4 contains 
the analysis of the column density and the temperature of 
the disk/ring of MMl, followed by a discussion in section 
5. Finally, in section 6 we draw our main conclusions. 

2. OBSERVATIONS 

The SMA observations were taken during two epochs: 
on 2005 August 24 in the compact configuration (giv- 
ing a synthesized beam of S'.'l x S'.'O) and on 2007 May 
27 in the very extended configuration (giving a syn- 
thesized beam of 0'.'7 x 0'.'4). During both epochs the 
receiver was tuned at 215/225 GHz. The phase cen- 
ter of the telescope was RA(J2000.0)= 18''19'"12U and 
DEC(J2000.0)= -20°47'30'.'0, and the correlator pro- 
vided a spectral resolution of about 1.14 km s^^ at the 
observed frequency (with 0.81 MHz of channel width). 
The co ntinuum data wer e edited and calibrated using 
Miriad (jSault et al.lll995| ). The flux uncertainty was es- 
timated to be ~ 20%. The detailed description of the 
observing setup as well as the calibration process of the 
continuum emission data can be found in Paper I. 

Nine emission lines were detected with the compact 
configuration (Table [Ij , four of which were also detected 
with the very extended configuration (SO2 I63 13-I62 14, 
SO2 132,12-131,13, SO 55-44, H2CO 3i,2-2i,i). The so- 
lutions of the self-calibration performed on the contin- 
uum data were transferred to the line data, and then the 
molecular lines were imaged, cleaned, restored, and ana- 
lyzed using the Miriad and AIPS (developed by NRAO) 
packages. The line data were also corrected for the half- 
channel error in the SMA velocity labeling discovered in 
2007 November. After that, we used a channel separa- 
tion of 1.136 km s~^ for all the lines. The calibration was 
performed in the same manner for the compact and very 
extended configuration data, except for the fiux calibra- 
tion step (see Paper I). Finally, the two data sets (com- 
pact and very extended configurations) were combined 
for all the lines giving maps with an intermediate angu- 
lar resolution (hereafter combined data). The average 
rms noise level is ~70 mJy beam"^ in the line chan- 
nel maps obtained with the compact configuration data, 
~85 mJy bcam^^ in those obtained with the very ex- 
tended configuration, and ^85 mJy beam-^ in those ob- 
tained with the combined data, all with a natural weight- 
ing. Throughout the paper, all the velocities are given 
in LSR. 

We use the high angular resolution data (very extended 
configuration) to analyze the disk kinematics and the 

^ The Submillimetcr Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 



low angular resolution data (compact configuration) to 
extract the spectra and estimate the disk temperature 
and column density. In addition, we obtain two differ- 
ent sets of images using the combined data to examine 
the extended emission with different angular resolutions. 
First, we set the robust parameter to 0.3. It yields a 
synthesized beam of about 1'.'2 x 0'.'9 (Table [1]). Second, 
we apply a taper of 182 kA and restrict the visibilities up 
to 200 kA, thus improving the sensitivity to the extended 
emission. The synthesized beam of the tapered data is 
about 4'.'0 X 2'.'1 (see Table H]). 

3. RESULTS 



3.1. C^^O 



emission 



Fig. 1 shows the integrated emission of the C-'^'^O 2- 
1 molecular line superposed on the first order moment. 
The C^^O emission extends mainly along a south-east 
to north-west lane of -- 23" x 10" (-- 0.2 pc x 0.1 pc) 
and a position angle of ^ 120°, approximately perpen- 
dicular to the thermal radio jet (position angle of 21°). 
The two main peaks of the C^^O 2-1 emission appear 
on both sides of the thermal radio jet (see Fig. [T|), with 
the strongest peak at a position between the two dom- 
inant (sub)milhmeter sources, MMl and MM2. MMl 
(the source from which the radio jet is launched) , is lo- 
cated at the bottom center of the C^^O 2-1 lane emission. 

Fig. [1] also shows that the vlsr of the line is almost 
constant along the lane structure. The line is narrow 
and has a symmetric profile, with a full width at half 
maximum (FWHM) of 2.8 km s'^ ±0.03. A Gaussian 
fit to this line profile gives also a central velocity of 
11.8±0.03 km s"^. In what follows we will consider this 
value as the systemic velocity. This velocity agrees well, 
to within our spect ral resolution, with that measured by 
I Gomez et all (|2003[) (12.2±0.1 km s^^) from the emission 
of the ammonia core. 

The critical density of the C^^O 2-1 line is relatively 
low, ricrit ~ 10'*cm-'^,so we expect that the emission de- 
tected by the SMA is likely thermalized. Assuming LTE 
conditions will yield a lower limit for the total mass of 
the extended molecul ar cloud. We assume a C^^O abun- 
dance of 4.7 X 10"* (jFrerking et al.|[T98^ and an cxci- 
tation tempe r ature of 24 K, equal to the T^ot found by 
I Gomez et al.l (|2003l ) using th e NH3 VLA obse rvations. 
From equations Bl and B2 of iFrau et al.l ()201Gf ) and for 
optically thin emission, we obtain an average molecular 
column density of 7 x 10^^ cm"^. Hence, we estimate a 
lower limit of 38 Mq for the total mass. 

3.2. 5*0 emission 

The SMA correlator setting includes two SO transi- 
tions: SO 78-77 and SO 55-44. We detect SO 78-77 only 
in the compact configuration due to a problem with the 
correlator at the frequency of this line in the very ex- 
tended configuration. The SO 78-77 line emission from 
the compact configuration data appears as an unresolved 
clump toward MMl. On the other hand, the emission 
from the SO 55-44 line at the systemic velocity of the 
dense core shows a clumpy distribution mainly to the 
north of MMl, having a conical shape which surrounds 
the north lobe of the thermal radio jet (Figs. [5] and [3]). 
Some of the weaker clumps seem to coincide with the 
direction of the north-west and south-east outflows ob- 



TABLE 2 
Gaussian pits to the low angular resolution data toward MMl. 



Line 



'-'peak 

K 



Vlsr 
km s~^ 



fwhm(^) 

km s^^ 



Gaussian area 
K km s^^ 



(b) 



K 



Zero- to-zero 
km s"-"^ 



S02 163,13-162,14 


1.32 ±0.03 


13.04 ± 0.05 


5.9 ±0.2 


8.2 ±0.4 


0.03 


5.5-19.0 


SO2 132 12~13l 13 


1.33 ±0.05 


12.9 ±0.1 


6.0 ±0.4 


8.5 ±0.4 


0.07 


6.5-18.0 


SO2 176,12-185.13 


0.56 ±0.05 


13.2 ±0.3 


5.7 ±0.4 


3.4 ±0.8 


0.06 


8.5-20.0 


SO2 202,18-193,17 


0.88 ±0.05 


12.8 ±0.2 


6.6 ±0.4 


6.2 ±0.7 


0.06 


6.0-19.0 


SO2 64,2-73,5 


0.49 ± 0.05 


13.4 ±0.2 


6 ±0.5 


3.3 ±0.7 


0.05 


6.5-20.0 


SO 55-44 


2.29 ±0.04 


12.86 ± 0.04 


5.65 ±0.05 


13.8 ±0.5 


0.04 


6.0-19.0 


SO 78-77('*) 


0.5 ±0.1 


9.8 ±0.3 


3 ±0.5 


1.3 ±0.4 


0.07 


6.5-18.0 




0.63 ±0.05 


13.8 ±0.3 


4 ±0.5 


2.5 ±0.4 






H2C0 3i,2-2i,i('i) 


3.6 ±0.1 


9.72 ±0.05 


2.1 ±0.1 


8.3 ±0.6 


0.10 


5.0-17.0 




3.7 ±0.1 


13.10 ±0.05 


2.6 ±0.2 


6.56 ±0.6 







Note. — The spectra were obtained over a box (-8,-8,7,7) ar csecs. Gaussian fits were carried 
out with a program based on the algorithm bv lCanto et al.ir2009l . 

(^) Calculated from cr„ ■ 2-\/ln2. 

f*^' It is obtained as ^/n ■ Speak • Uv ■ 

f'^' Estimated from the free-line part of each spectrum. 

t'*) A simultaneous two-gaussian fit was performed for this transition. 
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Fig. 7. — Image of the VLA radio continuum emission at 3.5 cm (red contours), SMA 1.4 mm continuum emission (black contours) 
and the SMA integrated emission of the SO2 132,i2-13i,i3 line (colour scale). The contour levels of the 3.5 cm map are 5, 6, 10, 15, 40 
and 200 times 8.8 /^mjy beam~^, the rms of this image. The contour levels of the 1.4 mm image are -5, 3, 3, 5, 10, 60 and 120 times 
0.003 mjy beam~^, the rms of this image. The SO2 132,i2-13i,i3 molecular emission is integrated over the velocity range 7-19 km s~^, 
and the wedge panel in the right shows the colour scale intervals. In the bottom left corner the synthesized beams of the cm image (red 
ellipse) and the mm images (black ellipse) are shown. 



served by iQiu fc Zhand (|2009[ ) . The redshifted emission 
(with respect to the cloud velocity) is compact and is 
mainly distributed around MMl and MC (Fig. [3]). The 
SO 65-44 line profile at the position of MC (top panel 
of Fig. 2]) has associated a clear redshifted wing with 
emission spreading up to a few tens of km s~^. 

3.3. H2 CO emission 

The emission of the H2CO 3i^2^2i,i line is more ex- 
tended than that of the SO 65-44 line, resembling the 
distribution of the elongated C^^O 2-1 emission around 
the systemic velocity (Fig. [5] and [6]). Furthermore, at 
the 11.9 and 13.1 km s~^ velocity channels the emis- 
sion appears to follow the direction of the two out- 
flows mentioned above. At redshifted velocities ( ulsr 
> 13.1 km s~^) the emission is mainly distributed around 
the north lobe of the radio jet, similarly to that observed 
in the SO redshifted emission. Contrary to the SO, the 
H2CO 3i^2-2i^i presents emission around MM 2. On the 
other hand and as occurs with the SO, the H2CO 3i^2^ 
2i_i emission is enhanced at the positions of MMl and 
MC. Toward MC, the H2CO 3i,2-2i,i line shows a high 
velocity red wing, more prominent than that observed in 
SO 65-44 line profile (bottom panel of Fig. |4]). 

3.4. 5*02 emission 

Five SO2 lines were detected with the SMA in its com- 
pact configuration (Table [2]). All of them show unre- 
solved emission at the origin of the thermal radio jet, 
which coincides with the MMl (sub)millimeter contin- 
uum source (Fig. [71). SO2 163.13-162,14 and SO2 132, 12- 
13i,i3 were also detected by the SMA in its very ex- 
tended configuration. The other SO2 lines were unde- 
tected due to their low intensity, their (probably) par- 
tially resolved nature and the lower sensitivity of the high 
angular resolution data. The panels a) and b) of Fig. [5] 
show the velocity channel cubes of the SO2 I63.13-I62.14 
and SO2 132,i2-13i,i3 lines (images of the combined data 
with robust 0.3). The molecular structure at the position 
of MMl appears at the 8.7 km s~^ to 16.6 km s""'^ velocity 
channels, while the peak velocity is vlsr — 13.0 km s~^ 
(see Table [2|) . The peak velocity of the molecular struc- 
ture is therefore redshifted with respect to the large scale 
dense core velocity by ^ 1.2 km s~^. The blueshifted 
channels (8.7 to 11.0 km s^^) are seen to the east, while 
the redshifted channels (14.4 and 16.6 km s~^) are seen 
to the west with respect to the thermal radio jet and the 
1.4 mm continuum peak position. 

The left panels of Fig. |9] show the zero order moment 
(i.e., integrated emission) and first order moment (i.e., 
integrated velocity weighted by the intensity) images 
built up from the SO2 I63.13-I62.14 and SO2 132,i2-13i.i3 
high angular resolution cubes. The molecular structure 
toward MMl is compact and its peak coincides with 
the position of the thermal radio jet within 0'.'2. The 
deconvolved sizes of the zero order moment images of 
SO2 163,13-162,14 and SO2 132,i2-13i,i3 are presented in 
Table m' Assuming a distance of 1.7 kpc to IRA S 18162- 
2048, the equivalent radius of the emission [^ Area / t:) 
is less than about 860 AU for both lines. The first order 
moment images clearly show the southeast-northwest ve- 
locity gradient nearly perpendicular to the radio jet axis. 
In fact, the position angle between the most redshifted 



and blueshifted channels are about ^ 110° and ^ 130° 
for the SO2 163,13-162,14 and SO2 132,i2-13i,i3 lines, re- 
spectively. In section §6.1 we argue that the SO2 emission 
arises from a disk/ring rotating structure. 

4. ANALYSIS 

4.1. RADEX 

In this section we estimate the physical characteristics 
of the MMl disk, such as volume density and tempera- 
ture by means of RADEX modeling. We used RADEX to 
simulate the line intensities of the five observed SO2 tran- 
sitions (Fig. [TU]). The RADEX code is a non-LTE molec- 
ular radiative transfer c ode which assumes an iso thermal 
homogeneous medium (jvan der Tak et aTl I2007D . This 
assumption is reasonable as a first approach to constrain 
the physical properties of the gas traced by the SO2 since 
the observed lines are optically thick (see below). We 
explored a range of values between 60 and 300 K in the 
kinetic temperature, between 10^ and 10^^ cm^'^ in the 
volume density, n(Il2), and between 10^** and 10^^ cm~^ 
in the SO2 column density, N{S02)- In order to con- 
strain the physical properties of the SO2 emitting region 
the next scheme was followed: 

• We first computed the x^ parameter (which is a 
measure of how well the model fits the observa- 
tions) for the line ratio of the 163,13-162,14, 176, 12- 
185,13, 202,18-193,17 and 64,2-73,5 SO2 transitions 
with respect to the 132,i2-13i,i3 SO2 transition. 
This approach assumes that all the transitions 
trace the same region (the images from the com- 
bined data of these lines give similar deconvolved 
sizes all with deconvolved major axis < 3"). Fig- 
ure [TT] shows in grey scale a set of the best so- 
lutions in the n(H2)-Tkin plane for four different 
values of the SO2 column density. The best set of 
solutions (darkest areas of Fig. [TT|) indicates that 
the SO2 arises from hot and very dense molecular 
gas, as expected if it is associated with the MMl 
disk. Indeed, the line ratios do not constrain well 
the physical values of the gas but rather yield lower 
limits: n(H2)> 3 x 10^ cm'^ and Tkin > 120 K. In 
addition, RADEX predicts a very large molecular 
column density of SO2, iV(S02)> 3xlOi'^cm~2, in- 
dicating that the observed lines are optically thick. 

• As an additional constraint, we take into account 
the possible values of the SO2 132,i2-13i,i3 and 
163,13-162,14 brightness temperatures. As a lower 
limit we can use the main beam brightness temper- 
ature measured in the highest angular resolution 
maps, Tmb = 60 K (it is shown in Figure [11] as a 
dashed line). However, it is not enough to further 
constraint the physical conditions of the gas asso- 
ciated with the SO2 emission. As a further con- 
strain, the true brightness temperature of the SO2 
can be estimated if the source size is known. The 
high angular resolution maps of the 132,i2-13i,i3 
and I63 13-I62 14 integrated emission yield a size for 
the SO2 of ~ iy.'S-l'.'O (see Table [3]). This size can 
also be independently obtained by comparing the 
main beam brightness temperature of the low and 
high angular resolution maps of these two lines. For 
a Gaussian distribution of the emission, the ratio 



TABLE 3 

Size and kinematics of the SOo disk. From the high angular resolution data. 



Line 



Deconvolved Size 



AU 



Vdisp 

km s 



(b) 
1 



©disp 

M0 



(c) 



Mdyn 



(d) 



SO2 163,13-162,14 
SO2 132,12-131,13 



0.88 ±0.08 
1.20 ±0.06 



0.81 ±0.05 
0.81 ±0.05 



130 ± 13 

165 ±6 



720 ± 50 
840 ± 50 



8.6 ±0.6 

6.7 ±0.6 



0.87 ±0.15 
1.03 ±0.15 



15 ±5 
11 ±3 



(^) Equivalent radius of the disk. For instance R — \Jhraai ■ 6,n,n/2. 

('=' Maximum velocity dispersion from the 50% contour of the p-v diagrams. The uncertainty is 

half the velocity channel. 

f'^' Maximum spatial offset from the 50% contour of the p-v diagrams. The uncertainty is half 

the minor synthesized beam axis, which i s roughly in the direction of the disk . 

("*) Obtained using the expression (see e.g. JPatel et al.ll2005l : [Beuther fc Walsh|[2008l : IZapata et a!] 
[2OOI '): (M/Mo)/sin^(i) = 1.121 x 10"^ x {{Av/km s"^)/2)^ x (R/AU), where i is the inclination 
of the disk with respect to the plane of the sky, Av is the velocity difference measured at the 50% 
contour level and R is the radius fitted with the AlPS task IMFIT. 
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Fig. 8. — SMA channel velocity images (combined data with robust 0.3) of several lines toward the position of MMl, the source from 
which originates the thermal radio jet. The lines are presented in the following order: row a) SO2 163,13-162,14, row b) SO2 132,i2-13i,i3, 
row c) SO 55-44 and row d) H2CO 3i,2-2i,i. The velocity channels are 1.14 km s~^ wide and the central velocity is indicated in the upper 
left corner of the channels. The disk velocity center is 13.0 km s~^ while the dense core velocity is 11.8 km s~^. Contours are -5, -3, 3, 5, 
7, 9, 11 and 13 times the rms noise level of each image, given in Table [T] The direction of the radio jet and the direction perpendicular 
to it are marked by two straight lines in each channel. The last channel of each row shows the zero order moment of each line (contours) 
superimposed on the 3.5 cm continuum radio jet image (grey scale). The synthesized beam sizes appear as a blue ellipse in the first channel 
of each row. It is evident that the peak of the emission moves from east to west on the SO2 velocity cubes. 



TABLE 4 
Size of the SO2 disk. IL From high to 

LOW ANGULAR RESOLUTION RATIO 



Parameter 


SO2 

132,12-13l,13 


SO2 
163,13-162,14 


tT (k) 

TiZ (K) 

^ mb / -* mb 

fi (arcsec) 

n (AU) 


50 ±8 


53 ±10 


1.33 ±0.05 


1.32 ±0.03 


37.6 ±6.2 
0'.'62 ± 0^10 
1054 ± 170 


40.2 ± 7.6 
0'.'59±0yil 
1003 ± 187 



of the low to high angular resolution line intensity 
plow _ ta2 , n2^//fl2 ,^2^^ ^^^^^ 



ihigh 



is T ^^^/T 



mb 



{91^ + n^)i{e\ 



ri is the source size, ^high and ^low are the beam 
sizes of the high and low angular resolution maps, 
respectively. Table |4] shows the intensities of the 
132,i2-13i,i3 and 163,13-162,14 lines for the low and 
high angular resolution maps. Using this method, 
the SO2 source size is slightly smaller, 0'.'6 ± O'.'l 
(Table H]). Therefore, we adopt a SO2 disk size 
in the range of O'.'5-l'.'O (425-850 AU). This im- 
plies that the brightness temperatures of the 182,12- 
13i,i3 and 163,13-162,14 lines are ~ 62-84 K and 
68-94 K, respectively. The area within the solid 
lines of Figure [Til shows the possible solutions that 
yield an SO2 132.i2-13i.i3 brightness temperature 
within this range. 
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Fig. 9. — SO2 163,13-162,14 and SO2 132,i2-13i,i3 zero (con- 
tours) and first order (colour scales) moment images (left panels) 
toward MMl. Their corresponding position- velocity diagrams with 
a perpendicular crosscut through the jet axis (P.A.= 111°) at the 
position of MMl are shown in the right-hand panels. The contour 
levels for the zeroth order moments are -30, 30, 50, 70 and 90% 
the flux peak, and the colour bar is in km s~^. The emission is 
integrated over the velocity range 7-19 km s^^ for both lines. The 
yellow dotted line marks the orientation of the radio jet. The yellow 
circle is centered at the position of the millimeter peak continuum 
emission of MMl and has the equivalent radius of the upper limit 
of the size obtained from the 1.4 mm millimeter image (Paper I). 
In the position-velocity diagrams (panels on the right) the vertical 
dotted line marks the disk velocity. East is at the top and west is 
at the bottom of both position-velocity diagrams. The contours 
are 50% to 95% with a step of 15% the peak fiux. 
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Velocity (km.s ) 

Fig. 10. — Spectra of several transitions toward MMl from low 
angular resolution data. The dashed line is placed at the velocity 
of the disk peak emission (13.0 km s~^ on average). Intensity is in 
Kelvin. 



By combining the two aforementioned criteria and us- 
ing the 99% confidence interval of the x^, the possible 
reasonable solutions get significantly reduced. Thus, for 
example, iV(S02) = 1x10^* cm~^ has better x^ solutions 
than lower values of the SO2 column density, but these 
solutions predict a brightness temperature significantly 
higher than the expected value. The best solutions are 
for a kinetic temperature of Tkin = 120-160 K and a vol- 
ume density of n(H2) si 2 x 10^ cm~'^. These values are, 
approximately, in agreement with those derived from the 
dust emission associated with MMl disk (Paper I). The 
possible SO2 column densities are in the 4-6 x 10^^ cm~^ 
range. 




Fig. 11.— Plot of the set of RADEX solutions in the n(H2)- 
Tkin plane for the emission associated with the MM 1 disk for four 
different values of the SO2 column density (from bottom-left to 
top-right panels): 3, 4, 6 and 10 X 10^ cm~ . The grey image 
gives the values of the x^ solutions for the SO2 line ratios (see 
S 14.111 . The increase of the grey tonalities (from lighter to darker) 
indicated a decrease in the x^ values. The solid and dashed white 
contours show the 68% and 99% confidence region of the x ■ The 
dashed black line shows the range of solutions for a brightness 
temperature of the 132,i2-13i,i3 line of 50 K, which is the value 
measured at the high angular resolution map. The area between 
the two black solid line shows the range of valid solutions for the 
132,i2-13i,i3 brightness temperature if the SO2 has an emitting 
size between 0'.'5 and I'.'O (see ^ 14.111 . 



The abundance of the SO2 emission can be estimated 
from the derived SO2 column density taking into account 
the filling factor of the observations. In Paper I, we used 
the 1.4 mm dust emission at an angular resolution of 
~ O'.'Sl to derive a beam averaged gas column density of 
iV(H2) ~ 1.1 x 10^^ cm~^. For this angular resolution 
the filling factor of a source with a O'.'5-l'.'O size is 0.49- 
0.80. Therefore, the SO2 beam averaged column density 
at this angular resolution is 2.0-4.8 x 10^^ cm~^, which 
implies a SO2 abundance of X[S02] ~ 1.7-4.2 x lO"*. 

5. DISCUSSION 

5.1. A rotating molecular disk/ring toward MMl 

5.1.1. Evidence for a disk 

The emission of the SO2, SO and II2CO at the position 
of MMl seems to consist of a compact rotating circum- 
stellar structure, probably a rotating gaseous disk or ring 
(Figs. OandlHl). AU the species detected with the SMA 
toward MMl are expected to be present in a ^ 150 K 
disk, in which a central prot ostar is e vaporat ing the gas 
from the grain mantles fe.g.. lCharnlev 1997; Maret et alj 
12004( 1. In addition, the SO2 abundance (2-4x10"^) ob- 
tained with RADEX is consistent with that found in the 
innermost parts of the envelopes of h ot cores, where the 
temp erature is greater than 100 K ()van der Tak et aTl 

\2om . 

The case for a rotating disk surrounding the proto- 
star in MMl becomes stronger in light of the new re- 
sults reported here. The velocity gradient within the 
molecular emission of the SO2 lines appears compact, 
with a radius R~425-850 AU, surrounding the compact 
(R< 150 AU) 1.4 mm dusty disk and the thermal radio 
jet (Fig. [7]). Furthermore, the position angle between the 
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extreme blue and redshifted channels of two SO2 lines 
(see §3.4), is roughly perpendicular to the radio jet axis, 
with a position angle 110-130°. 

5.1.2. Dynamical mass 

We now estimate the dynamical mass with the 50% 
contour of the SO2 163^13-162,14 and SO2 132,i2-13i,i3 
position- velocity diagrams in Fig. El We assume an equi- 
librium between the centrifugal and gravitational forces 
and also that the disk is seen edge-on. We use the 
virial equation, Mdyn = {v'^R)/G, which requires that 
the dominant force is due to gravity. This gives a dy- 
namical mass of 15±5 M© and 11±3 Mg (Table |3l) by 
using the SO2 163,13-162,14 and the SO2 132,i2-13i,i3 
data. The virial approximation is valid when the system 
is spherically symmetric regardless of its density distri- 
bution. It is in addition a good approximation for thick 
disks (as appears to be the case of MMl, §5.1.3). 

The continuum observations toward MMl (Paper I) 
have shown a very compact disk-like structure {R < 
150 AU, 4.1 M0) and a more extended and weaker enve- 
lope (i? ~ 1500 AU, 1.5 Mq). On the other hand the SO2 
rotating structure has a radius of 425-850 AU. For this 
radius range the gravitational effect of the compact disk- 
like structure could be considered as part of the central 
mass, therefore, using the virial equation is a reasonable 
approximation to find out the protostar plus inner-disk 
mass. 

To probe the impact of the extended molecular struc- 
ture on the centrifugal balance of forces in MMl, we 
consider an additional term due to a flattened disk, 
representing the SO2 disk-like structure. We assume 
that this structure has a mass of 1.5 M© and a ra- 
dius of 1500 AU, the characteristics of the dusty enve- 
lope of MMl. We also assume a disk surface density 
inversely proportional to the distance from the proto- 
star {p{R) oc R~^). In addition, we include another 
term to the balance of forces due to a pres sure gradi- 
ent {llp){dP/dR) (see iTorrelles et all [ToM) . To esti- 
mate this pressure term we use either a sound speed 
of 0.7 km s^^ (corresponding to a 150 K gas temper- 
ature), or a turbulent plus thermal random velocity of 
~ 1 km s^^. Both terms (extended disk and pressure 
gradient) could induce a departure from the Keplerian 
behavior of the disk. However, we find that both are 
one order of magnitude smaller than that of the central 
mass. Thereby, the central mass obtained is essentially 
the same as that estimated with the virial approxima- 
tion. 

If the SO2 rotating structure is not edge-on (i.e., 
iy^ 90°), the derived total mass is just a lower limit. Nev- 
ertheless, the edge-on approximation seems adequate for 
the MMl disk, as the thermal radio jet is one of the 
largest jets (~ 5.3 pc in projected distance) observed up 
to now. T herefore, it is expe cted to be close to the plane 
of the sky ([Marti et al.lll993l ). However. I Yamashita et all 
()1987() estimated an inclination between 60° and 74° by 
modeling the geometry of the infrared refiection nebula, 
which would increase the dynamical mass to a range of 
14-20 Mq. 

5.1.3. Stability of the disk 

Due to the high mass of the compact dusty disk (about 
4 M0), it is reasonable to ask whether the molecular 



disk is stable against gravitational disturbances. If we 
take as the disk parameters T=: 120-160 K, N= 1 x 
10^^ cm~^ (and hence the surface density, E = iV ■ /i • 
mn = 40 g cm~^, using a mean molecular weight of 2.3 
for molecular hydrogen and helium at 10 K), and the disk 
velocity 3.4-4.3 km s'^ at 876-740 AU (Table |31) from 
the center of the disk, we can approx imate the value of 
the Toomre parameter ()Toomrelll964[ ). defined as 

Q 



TT-G-T, 



In the equation above, Cs is the sound velocity, k is 
the epicyclic frequency (which can be approximated as 
n for Keplerian motion or as 2r2 for solid body rota- 
tion, n being the angular velocity), G is the gravita- 
tional constant and E is the surface density. If Q is 
less than 1 then the disk is unstable. First, we estimate 
the Toomre parameter using the thin disk assumption 
(e.g.. lKratter fc Matzner|[2006D . and a high epicyclic fre- 
quency (that of solid body rotation) . The resulting value 
is Q — 0.42-0.72. It appears that the molecular disk 
would be unstable even in the unlikely case of a rigid 

body rotation. 

Ho wever. iCesaroni et al.l (|2007l ) (see also lDurisen et al.l 
120011 ) pointed out that the thickness of disks should have 
a great stabilizing effect on these. Following their ap- 
proximation for thick disks by assuming a disk thickness 
H= c2/(7rGE) = 35-47 AU, a disk mass of Md = 4 Mq, 
and a total mass in a wide range of Mtot = 10-20 Mq, we 

obtain Q= y/2 ■ (H/R) -(Mtot/Md) = 0.40-0.80, which 
is again smaller than unity. 

Although the Q= 0.8 case is closer to unity, the MMl 
disk-like structure seems therefore probably unstable. 
In such a case, the gravitational perturbations can lead 
to the appearance of dense spiral arms, rings, arcuate 
structures or even dense fragments, which could result in 
the fo rmation of a clumpy structure (iDurisen et al.|[200ll 
[2OOI . 

5.1.4. Mass of the protostar 

We can also roughly estimate the mass of the proto- 
star at the position of MMl. If we assume an edge-on 
disk (which provides again with a lower limit), since the 
gas plus dust mass of the MMl disk was estimated from 
the 1.4 mm millimeter continuum emission (see Paper I) 
to be about 4±1 Mq, the mass of the protostar would 
be about 7±4-ll±6 Mq. Although the estimation has 
an uncertainty larger than a 50%, this agrees well with 
the expected mass for a > 3300 Lq (Paper I) protostar, 
whi ch could be related to a -- 10 Mq ZAMS star (Table 
5 of lMolinari eTaI1fT998[) . 

5.1.5. Absorption of disk emission by a foreground cloud 

The blueshifted emission of the SO 55-44 and 
H2CO 3i,2-2i.i lines from the disk appears to be missing 
(Fig. [8|). One possibility is that the interferometer has 
filtered out some extended fiux. However one would ex- 
pect some remnant emission from the rotating structure. 
Another possibility is that gas from the foreground com- 
ponent of the cloud between 8.5 and 11.8 km s~^ (Figs. [3] 
and in]) may have absorbed the blueshifted emission from 
the disk. Firstly, the extended dense molecular compo- 
nent traced by the C^^O and H2CO appears blueshifted 
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at the same velocity range where the missing velocity disk 
component is observed in SO 55-44 and H2CO 3i,2-2i^i. 
And secondly, the other SO transition, SO Tg-Ty, shows 
only compact emission toward MMl, as in the other SO2 
lines. This SO line and also the SO2 lines, have a higher 
excitation, therefore they can only be excited in the hot 
environment around the disk. 

5.2. Emission from the cavity walls 

Figs. [5] and [3] show the SO emission extending about 
10"-15" to the north of the position of MMl. The SO 
emission appears clumpy, with some of the clumps proba- 
bly associated with the molecular outflows of the region. 
In particular, the SO 55-44 line emission, and also in 
a lesser measure the H2CO, seems also to surround the 
thermal radio jet axis, with a geometry reminiscent of 
a V-shape. A possibility is that part of the emission of 
these lines is tracing the walls of the outflow cavity exca- 
vated by the thermal radio jet originated at the position 
of MMl (as that of outfl ow caviti es associated with low- 
mass protostaJS^_see^e.g. , |L^^^L2000;^Ajx^^ 
l2OO6Hj0rgensen et a l. 2007:' Santiago-Garcia et al.ll2009D . 
This quiescent molecular emission (i.e., showing narrow 
lines) associated with powerful outflows could be photoil- 
luminated by the UV radiation generated in the shocks 
within t he outflows, as has been observed in several HH 
objects (|Girart et al.l[200l [200l[Viti et al.l[200l[2006h . 
Indeed, in the case of the HH 80-81 and HH SON system, 
there is a shock-induced photodissociati on region along 
all th e thermal radio jet, reported by iMolinari et al.l 

(pool . 

5.3. Is the Molecular Core actually a hot core? 

The only molecular lines clearly detected toward the 
molecular core (MC) northwest of MM2 are SO 55-44 
and H2CO 3i^2"2i^i. This source has been previousl y 
identified as a hot molecular core (|Oiu fc Zhanel [20091 ) . 
Both lines show a narrow component (e.g., the FWHM 
of the SO line is ~ 2.5 km s^^) centered at wlsr —10- 
12 km s~^, and a redshifted wing spreading up to about 
30 km s^^. The narrower component is related to the 
quiescent gas in the molecular core, while the component 
with the wider velocity range is probably due to shocked 
gas entrained by the outflows of the region. Herein, we 
discuss the possible origin of the MC. 

It has been proposed that the MC is a hot core sur- 
rounding a massive protostar. This is based on the de- 
tection of hot core molecular tracers, s uch as CH3OH, 
HNCO, and OCS (|Qiu fc Zhand [2009h . However, we 
note th at the temperature estimated in the MC is ~ 45 K 
(jOiu fc Zhang 2009:) which is about a factor of 2 lower 
than the temperature found in hot cores. In addition, the 
lack of dust emission associated with the MC (< 10 mJy 
at 1.4 mm) implies that for the aforementioned temper- 
ature the tot al mass of the MC has an upper limit of 
~ 0.25 Mq (jQiufc Zhang|l2009[ ). too smaU for a hot 
core. Could it be pos sible that MC is a hot corino such 
as IRAS 16293-2422 (|Rao et al.ll2"009l )? At the distance 
of IRAS 18162-2048, IRAS 16293-2422 source B would 
have a flux density of 10 mJy at 1.4 mm, below the upper 
limit of our observations. But, the MC temperature is 
still lower than that of this hot corino and the intensity of 
the observed molecular lines in source IRAS 16293-2422 



B would be far from being detected with the current sen- 
sitivity at the distance of IRAS 18162-2048. 

A possible alternative is that MC is tracing the interac- 
tion of the dense molecular core with a strong shock pro- 
duced in a molecular outflow in a way that the external 
heating and the r adiation of the gas produced a strong 
outgassing phase (jViti et al.ll2006[) , exciting many molec- 
ular transitions, as it has been observed in other sou rces, 
suc h as Orion (iLiu et al.l 120021 : iZapata et al.l l2010t see 
also lChernin fc Wrightlll99fiir ^ 

There is a class I CH3OH maser (jKurtz et all 120041 : 
lOiu fc Zhanil2009f ) detected about 1" south-east of MC. 
This cla ss of masers is ass ociated with shocked molecu- 
lar gas (jKurtz et al.ll200^ . The most likely scenario is 
that the molecular outflow from MM2 has encountered a 
molecular core at the MC position. The northwest lobe 
of this outflow is redshifted and its path is crossed by the 
MC (see Figs. [3] and |6]). Therefore, it seems possible that 
the outflow could partially hit the molecular core at MC, 
producing a stationary bow shock which could entrain 
the molecular gas, as found in some molecular outflows 
such as L1157 and IR AS 04166+2706 (Zhang ctalJIiMl; 
iSantiago-Garcia et al . 2009; Tafalla et al. 20 ij). Given 
the location of the MC, it is also possible that it is not 
an isolated core, but it forms part of the cavity walls 
associated with the thermal radio jet (section §5.2). 

6. CONCLUSIONS 

We have performed an SMA observational study of the 
molecular gas toward the central region of IRAS 18162- 
2048, and have obtained the following main results. 

• The SO2 163^13-162,14 and SO2 132,i2-13i,i3 data 
show a compact (R~ 425-850 AU) molecular disk- 
like structure spatially coincident with the position 
of the compact millimeter dusty source MMl and 
the radio jet located at the IRAS 18162-2048 posi- 
tion. The molecular structure is perpendicular to 
the radio jet axis and it shows a clear velocity gra- 
dient that we interpret as rotation. The dynamical 
mass (11-15 Mq) allowed us to derive the mass of 
the protostar (7-11 Mq), making use of the mass 
that was inferred from the continuum data (Paper 
!)• 

The RADEX analysis constrains the disk-like 
structure temperature between 120 and 160 K and 
its volume density as > 2 x 10^ cm~'^. We have also 
found that the rotating system could be unstable 
due to gravitational disturbances. 

• The C^^O 2-1 emission shows an elongated struc- 
ture of about 0.2 pc x 0.1 pc with an orientation 
roughly perpendicular to the radio jet. This emis- 
sion seems to belong to the dense core in which the 
protostars associated with the MMl, MM2 and MC 
sources are forming. 

• The H2CO and the SO emission appears clumpy 
and is mainly enhanced toward the position of 
MMl and MC. It is possibly tracing in part the 
cavity walls excavated by the thermal radio jet. At 
the position of MC, it shows a line profile with a 
red wing extending up to tens of km s~^. This line 
profile is probably due to shocked gas entrained by 
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the outflows in the region. We speculate that a pos- 
sible scenario is the interaction between the outflow 
associated with MM2 and the dense molecular core 
at this position. 
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